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ABSTRACT 

We present the results of attempts to detect the ehipticity of dark matter halos using 
galaxy-galaxy weak lensing with SDSS data. We use 2 020 256 galaxies brighter than 
r = 19 with photometric redshifts (divided into colour and luminosity subsamples) as 
lenses and 31 697 869 source galaxies. We search for and identify several signal contam- 
inants, which if not removed lead to a spurious detection. These include systematic 
shear that leads to a slight spurious alignment of lens and source ellipticities, intrinsic 
alignments (due to contamination of the source sample by physically-associated lens 
source pairs), and anisotropic magnification bias. We develop methods that allow us 
to remove these contaminants to the signal. We split the analysis into blue (spiral) 
and red (elliptical) galaxies. Assuming Gaussian errors as in previous work and a 
power-law profile, we find fh — eu/eg = 0.1 ± 0.06 for red galaxies and —0.8 ± 0.4 for 
blue galaxies using transverse separations of 20-300 /i~^kpc, averaged over luminos- 
ity. Inclusion of the more realistic non-Gaussian error distributions and of the NFW 
density profile (which predicts much smaller ellipticity of the shear for scales above 
the scale radius) yields 0.60 ±0.38 for ellipticals and —1.4^2 for spirals. While there 
is no concrete detection of alignment in either case, there is a suggestion in the data of 
a positive alignment in the brightest lens sample of ellipticals. Our results appear to 
be mildly inconsistent with a previously reported detection by Hoekstra et al. (2004) , 
but more data and further tests are needed to clarify whether the discrepancy is real 
or a consequence of differences in the lens galaxy samples used and analysis methods. 

Key viTords: galaxies: haloes - gravitational lensing. 



1 INTRODUCTION 

Dark matter halo ellipticity, a robust prediction of ACDM 
according to N-body and hydrodynamic simulations, can 
in principle be detected using galaxy-galaxy weak lensing. 
In this paper, we attempt to detect the projected elliptic- 
ities of the dark matter halos of lens galaxies using data 
from the Sloan Digital Sky Survey (SDSS). Besides impli- 
cations of such a detection for CDM and theories of hierar- 
chical structure, a positive detection of halo ellipticity may 
also be used to constrain Modified Newtonian Dynamics, or 
MOND CMile:rom..l98a : .Sanders..l986.: Sanders fc McGauehl 
|20o2), an alternative theory of gravity that describes ro- 
tation curves well without dark matter, but predicts an 
isotropic lensing signal around isolated galaxies at scales on 
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which there is no baryonic matter (though a combination of 
MOND and dark matter cannot be ruled out by this test). 

Galaxy-galaxy weak lensing is a useful tool for studying 
the structure of dark matter halos on large scales. Because 
it is sensitive to all matter in the lenses, not just baryons, 
the signal (averaged over many lens and source galaxies) 
can be detected well beyond the extent of the light profiles, 
without the need for physically associated tracers such as 



been well-detected 


in several surveys (Brainerd et al.l 


199d; 


Hudson et alJ 


199J 


; 


Fischer et al.ll200a iSmith et all 


200 J; 


McKav et alJ 


2001 


iGuzik & SeliakI 120021: iHoekstra et al 


2004 


Sheldon et alj|2004l: iHoekstra et al.ll2004l:ISeliak et al 


|2005a 


), with increasing statistical precision and gradu- 



liooii), and therefore may also be useful as a tool to study 
dark matter halo ellipticity. 
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EUipticity of dark matter halo profiles has 
been predicted in CDM N-body simulations (e.g., 
iDubinski fc CarlberdlTool . and observed w ith non-lensin g 
methods on scales < 20 kpc (for a review, see ISackettll99sl) . 
Simulations typically predict triaxial halos with b/a > 0.7 
and c/a ^ 0.5, where b/a is the ratio of intermediate to 
long axis, and c/a is the ratio of short to long axis. The 
resulting mean projected ellipticity is then 0.3. It should be 
noted, however, that simulations that include baryons show 
reduced halo ellipticity in the inner regions, w ith A(c/a) 
and A{b/a) ~ 0.1-0.4 jKazantzidis et ahllioM) : the effect 
of baryons decreases with radius. One measurement of halo 
ellipticity usi ng galaxy-galaxy wea k lensing has already 
been reported. iHoekstra et al.l ^2004^ argued for a detection 
of ellipticity by measuring fh = Ch/eg, the ratio of the dark 
matter halo ellipticity to the ellipticity of the light profile in 



projection, and found fh ~ 0.77" 



+0.18 



nominally ruling out 



an isotropic lensing signal at the 99.5 per cent confidence 
level. 

A measurement of dark matter halo ellipticity naturally 
depends on alignment between the elliptici ty of the light 
distribution and of the matter distribution. iNavarro et alJ 
i2Q04l) find, using gasdynamical simulations, that for spiral 
galaxies on small scales (within the virial radius) the disk 
spin axis is aligned with the minor axis of the dark mat- 
ter halo, so when viewing the disks in projection we should 
see positive f^. On larger scales (1-2 h~^Mpc), when con- 
sidering the alignment of the ellipticity of the light of spiral 
galaxies relative to large-scale structure, the reverse appears 
to be true, so that fh should be negative on those scales. 
In contrast, comparison of N-body simulations of satellite 
distributions wit h the observed satellite distribution around 
the Milky Way llLibeskind et al.ll2005l IZentner et ai]|2005ll 
suggest that its disk spin axis is aligned with the halo major 
axis, thus predicting a negative fh on all scales. In short, 
predictions for spiral galaxies are conflicting. For elliptical 
galaxies, the halo and light ellipticities are expected to be 
aligned within the virial radius in the simplest models, giv- 
ing positive fh ■ 

There are, however, a number of difficulties inherent in 
these measurements. Some are theoretical; for example, if 
the dark matter and light ellipticities are not well aligned, 
then any tendency for halo ellipticity will be undetectable. 
Likewise, as will be shown in t|2.3l for some reasonable halo 
profiles, the ellipticity of the weak lensing signal is much 
less than the ellipticity of the projected matter distribu- 
tion, so the signature of halo ellipticity in weak lensing is 
quite small. However, several effects are under our control: 
the statistics (i.e., finding a large enough sample that halo 
ellipticity effects, will be statistically significant); and PSF- 
related systematics, which may contaminate the halo ellip- 
ticity measurement, can be studied and understood. 

The SDSS provides an excellent dataset on which to 
carry out these measurements by the above criteria. With 
3 X 10^ nearby {z < 0.25) galaxies with spectroscopy that 
can serve as lenses, or eight times as many bright (r < 19) 
lenses if we do not require spectroscopy, and roughly 3 x 10^ 
fainter galaxies to serve as sources, the SDSS should provide 
the statistical power to measure these effects. Also, as will 
be described in H3.3I we have the tools to isolate the effects 
of PSF systematics on the halo ellipticity measurement. 

We begin by outlining the weak lensing formalism and 



the models used for the ellipticities of the dark matter halo 
in Next, in Sj^l we describe the data used for these mea- 
surements and systematic effects that may bias the mea- 
surements. ^ outlines the method of analysis that we use 
to relate results from the data to our models in ^ Results 
are presented in ^ and the implications and suggestions for 
future work are in 

Here we note the cosmological model and units adopted 
for this work. Pair separations are measured in transverse 
comoving h~^kpc using the angular diameter distance (in a 
flat ACDM cosmology with Qm = 0.3, Qa = 0.7) at the lens 
redshift. In the units used, our results are independent of 
Ho. 



2 THEORY 

2.1 Galaxy-galaxy weak lensing formalism 

In this paper, we use the notat ion and approach to t he lens - 
ing signal computation from iMandelbaum et al.l (12005311 . 
hereinafter M05. In this section we briefly mention the main 
points of the analysis. 

Galaxy-galaxy weak gravitational lensing is sensitive to 
the difl^erential surface density AE(r) = E(< r) — E(r-) (for 
an axisymmetric halo) averaged over many lenses. This sur- 
face density can be related to a product of a tangential shear 
term and a redshift-dependent term, AE = 7tEc, where the 
critical surface density is defined in comoving coordinates 
via 



A-kG Di,Di{l + zi) 



(1) 



where D is the angular diameter distance, and / and s refer 
to the lens and source, respectively. 

The tangential shear is computed by finding the average 
value of tangential ellipticity, but the two are related by a 
factor of 2 according to our definitions; furthermore, we must 
divide by a shear responsivity factor. We write 



AE = 



Wi{etT.c)i 
~2^^;^ 



(2) 



where the summation is over lens-source pairs, TZ is the shear 
responsivity which des cribes the effects of a she ar on the 
tangential ellipticities jBernstein fc Jarvi j l20o3l . and the 
weights are defined by 



(3) 



1 9 ' 



where the weighting is by inverse measurement noise plus 
shape noise. Since typical rms ellipticities are of order 0.36 
and typical shears are of order 10"'', we must average over 
millions of lens-source pairs to get reasonable signal to noise. 
The shear responsivity and weighting scheme are discussed 
in more detail in M05. Once these sums have been accu- 
mulated, we must subtract off the signal measured around 
random points, and boost by the number of pairs relative to 
the number around random points to account for dilution of 
the signal by physically-associated pairs. 
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2.2 Our Model 

Here we describe the model used in this paper 
for halo ellipticity. This f ormalism was developed by 
iNataraian &i RefregieJ (I2OO0I) and applied to the elliptical 
isothermal sphere density profile. 

We attempt to determine a parameter / which is related 
to the ellipticity of AE, the differential dark matter surface 
density of galaxies. We model AE as having an isotropic 
component AEiso and some azimuthal variation as follows: 



AE™,dei (r) = AE.so (r) [1 + 2/69 cos (2 A^)] 



(4) 



where Eg is the observed ellipticity of the light distribution 
of the lens and A6' is the angle from the lens major axis. We 
can determine AEiso and / by minimizing 



X 



™i(AE„ 



AEO 



(5) 



where the summation is over lens-source pairs denoted by i 
with weights described in H2.1l and AEi = etEc for the pair 
i. By minimizing with respect to / and AEiso, we obtain 
the following joint solution for the two quantities^: 



AEiso (r) = 
/AE,so(r) = 



J2i TOjAEi 

J2i WiAEieg,i cos (2A0i 



(6) 



The solution for AEiso (r) is thus as expected, the usual aver- 
aging of tangential ellipticity over all pairs, and the calcula- 
tion of /AEiso (r) is a summation weighted by Cg cos (2A6'), 
which will pick out any tendency for AE to be larger 
(smaller) along the major axis of the projected light distri- 
bution, yielding a positive (negative) value for /AEiso (?"). 
This procedure can be used even if / is not a constant with 
radius, allowing us to measure f{r)AT.iso{r). 

We also will find it necessary to calculate AE45, the sig- 
nal with source ellipticities rotated by 7r/4. While this quan- 
tity must average out to zero by symmetry when averaged 
over A9, we will see in t|2.3l that it can have a sin (2A6') de- 
pendence. If we also rotate the lens ellipticities in the same 
direction as the source ellipticities, that sin {2Ad) depen- 
dence becomes a cos2A0 dependence as for the tangential 
shear, i.e. 

AE45,w.i(r) = AE,so(r) [2/4563 cos (2(A6l + 7r/4))] (7) 
so we can write 

r r . E» t»» AE..4569.. cos [2(Ae. + 7r/4)] 

/45AE,so(r)= 2E.«^.e^,.cosM2(Ae.-f./4)] ' 

Once these sums have been calculated, the usual pre- 
scription must be followed as already described: divide by 
2TZ, subtract off the random catalog signal from AEiso, and 
multiply by the boost factor B{r), the ratio of lens-source 
pairs using real lenses to the number using random points 
as lenses, which goes to one in the limit of no physically as- 
sociated sources included in the sample, and is larger than 



^ This results in minimization if uiiEg i cos(2A6i) is negligi- 
bly small, as would be expected if the sources are isotropically 
distributed relative to the lens major axis. The validity of this 
approximation is extensively discussed in i|3.c{.2l 



that at small transverse separations. Multiplying by this fac- 
tor corrects the signal for dilution due to these non-lensed 
galaxies included in the source sample. 

There are a number of difficulties to consider when at- 
tempting to detect dark matter halo ellipticity that compli- 
cate this simple model we have presented here, in particular 
difficulties in getting adequate signal to noise, the question of 
whether the light and mass ellipticities are actually aligned, 
choosing radial ranges to best carry out this measurement, 
and systematics that obscure the measurement. All of these 
issues will be addressed in the sections that follow. 



2.3 Connection to theory and other work 

In order to facilitate a comparison of our results with theor y 
and with observational results from iHoekstra et al.l i2004l) . 
we must relate the quantities we measure to theoretical pre- 
dictions for various halo density profiles. 

We consider several different types of profiles. The sim- 
plest is the general power-law profile, for which we derive an- 
alytic expressions for / correct to first order in the ellipticity. 
More realistic profiles, particularly the Truncated Isother- 
mal Sphere (TIS) and the Navarro- Frenk- White (NFW) pro- 
file, are also considered, with numerical results given. In all 
cases, we introduce ellipticity by replacing r in the isotropic 
form for Hi{r) with an elliptical coordinate to compute (an- 
alytically or numerically) the potential and, ultimately, the 
shears 7t(r) and 745 (r). We find that the degree to which 
the ellipticity of n is expressed as ellipticity of AE depends 
on the shape of the profile. 

For a dark matter halo with power-law space density 
profile p(r) oc r"*-^"*""^ and ellipticity independent of radius 
6h, oriented with major axis along 6 = Q and tt for simplicity, 
we can write the surface density as 



K{r,e) ^A\r{l - ycos (261) 



= Ar 



[1 + ^ cos {26) 



(9) 
(10) 



This definition ensures that for major axis a and minor axis 
fe, the relationship ch = {a^ — h'^)/{a? + b'^) holds. The next 
order term, proportional to cos^ (26), has a coefficient that 
is smaller than that of the cos (2^) term by a factor of 6h(l-f 
q)/4, which for typical values a ~ 1 and eh ~ 0.3 is roughly 
0.15. Consequently we will ignore all terms of order higher 
than cos 29 for the remainder of this calculation. 

For an axisymmetric (around the line of sight) mat- 
ter distribution, we could simply obtain AE by finding 
the average value of n within radius r and using AE = 
Ec [k(< r) — ^(r)]. However, for non-axisymmetric halos, we 
must do the full calculation, first determining the lensing 
potential corresponding to the surface density via 



V<I>=i 
2 2 



dr 



9$ 
dr 



(11) 



With an ansatz of $ = Em <I>o(m)r cos [mO), our first 
order approximation (Eq.|UJ gives a potential of 



2A 



{a - 2)2 



l+ ^^'-y cos(2g) 
a{a - 4) 



(12) 



where B — ehce/2. Since this expression diverges for a = 2 
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{p tx r~'^), it cannot be used to derive expressions for the 
shear and / for NFW profiles well beyond the scale radius. 
We can then use the potential to determine 7t and 745 

via 



1 



+ r 



dr 



d9^ 



(13) 



-1 d'^ ^ _ 

-r h r 

dedr 



,9$ 



For a general power-law potential, we therefore find 



It 



745 



-Ar 



1 + 



(Q-2)(a2-2a + 4)B 



-Ar 



Qf2(a — 4) 
4(2 - a)(l - a)B 



cos (261) 



(14) 



a^{a - 4) 



sin (261) 



We can then make the identification A'Eisoir) = 
T,caAr~°' /{2 — a) to relate back to our model in Eq. |1] 

In the special case of an elliptical SIS (a = 1), we have 
d^^/dr^ = 0, 7t 



K, and therefore 



7t = Ar 
745 = 



1 + y cos (261) 



(15) 



Comparison against our model in Eq.0] shows that our / — 
efi/4eg, so that for — eg {fh ~ 1) we expect / = 0.25. 
Dividing the measured / by the predicted SIS value f /fh = 
0.25 then allows us to compute the measured fh from our / 
value assuming the SIS density profile. 

More generally, for a power-law non-SIS model, we find 



f _ fh 

•' 4 



[a - 2)(a^ - 2a + 4) 



a{a — 4) 



(16) 



so for a = 0.8-0.9 (a typical value in the actual data at 
small transverse separations) / ~ 0.33/h. 

We are also interested in computing /45 for these pro- 
files. By comparison with previous equations, we find 



fis = fh 



(2-a)(l-a) 



>:{a — 4) 



4(q - 1) 
-2a + 4 



f (17) 



Consequently the expected signal with lenses and sources 
rotated by 7r/4 can be predicted relative to the unrotated 
signal, and the comparison of the two is an important con- 
sistency check. For a ~ 0.8-0.9 as is found in the data, 
f^hlf ~ —0.26 to —0.13, i.e. the rotated signal is smaller 
than the unrotated one and with opposite sign. As will be 
described in H3.3I to eliminate contributions from systemat- 
ics to /AE, we actually will measure / — /45, which for the 
general power-law profile is 



_ (q-2)(q-4) _ fh 



(a - 2)2 



(18) 



We also consider non-power law density profiles. 
The first such profi le that we will consider is the TIS 
ijBrainerd et aljfl996t l. which takes the form 



p(r) oc 



r2(r| -)- r^) 



(19) 



(note that what we call here is often called s; we use 
for simplicity of notation since it also appears in the NFW 



profile). This profile yields k oc for r <^ and oc r"'' 
for r 2> Ts. While analytic expressions can be derived for 
the shear in the spherical case, numerical integration must 
be used when r is replaced by an elliptical coordinate. For- 
tunately, K and 7 can b e computed for this model using 
GRAVLENS^ iKeetonll200ll) via subtraction of the shear from 



model 



an elliptical k oc (r^ + s^)-^/^ model from a 
(due to linearity of the shear), with r being replaced by an 
elliptical coordinate. We can use the predictions for those 
quantities to compute / — /45 by finding, for a profile with 
ellipticity e, 

A(r) EE 7t(r, e = 0)~'yt{r,e = n/2) 
A45(r) = 745(r,6» = 7r/4) - 745 (r, 6^ = 37r/4) 
T(r) = 7t(r, 6 = 0) + 'yt{r,e = n/2) 
A(r) 



(20) 



fir) 



fMr) = 



2eT{r) 
A45(r) 



2eT(r) 

The plot in Fig.^shows the results for /, /45, and /— /45 
as a function oir/vs for the TIS over a wide range of scales. 
(This plot was computed using e = 0.02, for which the vari- 
ation of the shear with azimuthal angle is dominated by the 
first-order cos {26) term; at more realistic values like e = 0.3, 
higher order terms may contribute as much as 18 per cent of 
the value of this first order term. However, as will be shown, 
since our results include statistical uncertainty larger than 
this value, we will henceforth neglect higher order terms.) 
Notice that / for the TIS is a declining function of trans- 
verse separation, and for the smallest scales shown on the 
plot approaches the SIS prediction f / fh = 0.25. Because /45 
is so large, / — /45 is quite close to zero for r ~ 0.6rs to the 
largest scales shown on the plot, r ~ Z rs- 

We also consider the NFW profile (|Navarro et alligS^ , 
which takes the form 



P = 



(r/r»)(l + r/r,)2^ 



(21) 



where rs = r^/c (concentration parameter). While the shear 
for t he spherical NFW mode l can be computed analyti- 
cally ijWright fc Brai nerd"200n'') , the additional complication 
of ellipticity necessitates the use of numerical integration to 
get the shear. So, we again use 7 from GRAVLENS to com- 
pute / and /45. A plot of f(r) for this model is shown in 
Fig. |5| As shown, the NFW model gives decreasing /(r) as 
for the TIS model. While at r = rs, the NFW model gives 
f/fh ~ 0.3 (slightly larger than the SIS) and /45 < /, / 
decreases at larger radii and /45 increases so that / — fiz, is 
quite small for r larger than about 2rs. 

For the sake of clarity. Table summarizes the defini- 
tions of all ellipticity-related parameters used in this paper. 



3 DATA 

The data used for thi s paper come fro m the Sloan Digital 
Sky Survey, or SDSS llYork et ahlbOOol) . an ongomg survey 



The latest version may be found via the link at 
\Utp://cfa- www. harvard . edu/ castles/ 
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f{r)/fi for TIS profile 




Table 1. Definitions of all measures of projected ellipticities (of 
dark matter halo and source distributions) used in this paper. 



/ 
fh 

/45 
f-f 

fint 



feff 



Ellipticity of AS = ASiso[l + 2/eg cos (2A6»)] 
Ellipticity of halo relative to light distribution, 

Ellipticity of AS45 = 2/45 ASisoEs sin (2A6I) 
Ellipticity of shear before boosting; 
7^ / if boost has azimuthal dependence 
Ellipticity of distribution of 
physically-associated "sources" 
Ellipticity of magnification bias 
5N/N oc K oc (1 + 2f„,ageg COS {2Ae)) 
Effective ellipticity of fi„t , fmag 
taking into account dilution by isotropic 
source distribution 



0.01 0.1 1 

r/r^ 

Figure 1. Plot of predicted f{r)/fi^ and f45{r)/fh for an ellip- 
tical TIS density profile dark matter halo; the maximum radius 
shown, 3rs, is larger than the scales used in this paper. Horizontal 
lines indicate the SIS predictions f/fh = 0.25 and /45//h = 0. 

f(r)/f^ for NFW profile 



M 1 1 1 




'{'■)/'h 






f„(r)/f, 




\ \ 


-[f(r)-f„(r)]/f. 








ml , 







0.1 1 10 100 

r/r^ 



Figure 2. Plot of predicted /{r)//^ and f45{r)/fh for an el- 
liptical NFW density profile dark matter halo. Horizontal lines 
indicate the SIS predictions f/fh = 0.25 and fif,/ fh = 0. 

that will eventually image approximately one quarter of the 
sky (10,000 square degrees). Imaging data is taken in drift- 
scan mode in 5 filters, uq riz, centred at 355, 469, 617, 748, 
and 8 93 nm respectively jFukugita at alJll996llSmith et alJ 
l2002t) using a wid e-field CCD JGunn et alJll998^ with pho- 
tometric monitor iHogg et al.||200 l|). After th e computation 
of an astrometric solution Jpier et al.ll2003l) . the imaging 
data are processed by a sequence of pipelines, collectively 
called Photo, that estimate the PSF and sky brightness, 
identify objects, and measure their properties. The soft- 
ware pipeline and photometric quality assessment is de- 



scribed in llvezic et a l.' 12004*). Bright galaxies and other in- 
teresting objects are se lected fo r spectroscopy according to 
specific criteria jEisenstein et alJi200ll : IStrauss et alJl20o3 : 
iRichards et alJl2002^ . The SPSS has had four major data 
relea ses: the Early Data Release or EDR (|Stoughtoi^^^l 
120021), DRl fAbazaiian ct al."2003'), DR2 llAbazaiian et al, 
2004), and DR3 ( Abazaiian ct al. 200^. 

The lens and source catalogs are very similar to those 
described in M05, except for the inclusion of the fainter 
lenses with photometric redshifts. Thus the description here 
will be brief, with the exception of the few differences from 
that work which will be described in more detail. 



3.1 Lens sample 

While M05 used only 3 x 10^ spectroscopic galaxies (r < 
17.77) as lenses, this work uses a larger sample of 2 million 
lenses as faint as r = 19 (model m agnitude) with photom et- 
ric redshifts from kphotoz v3_2 iSlanton et alJl2003al) in 
the range 0.02 < z < 0.5. There are several reasons for this 
change. First, in order to detect a possibly small azimuthal 
variation in the tangential ellipticity, we need a much larger 
sample of lenses. Because of the higher number density of 
lenses when we include these fainter lenses, and because of 
the larger photometric area coverage, the use of a fainter 
flux limit increases the size of our lens sample by a fac- 
tor of roughly eight. Second, while a lensing analysis aimed 
at determining the signal amplitude suffers from significant 
calibration uncertainty whe n galaxies without spectros copic 
redshifts are used as lenses iKleinheinrich et alJl2005l) . this 
work is primarily concerned with the ratio of /AS to AE, 
so errors in signal calibration are irrelevant. 

There is one possible error in the detection of / that 
could be affected by the use of photometric redshifts. The 
photometric redshifts are used for each lens to determine 
Da{z) and consequently the transverse separation between 
lenses and sources. If / is a constant value independent of ra- 
dius for the projected mass distribution, then errors in the 
computed transverse separation do not matter, since they 
will affect AE and /AE in the same way. However, if / 
decreases with radius, and errors in the photometric red- 
shifts tend to go in one direction (nonzero average bias), 
then the measured / can be systematically affected. If the 
photometric redshifts tend to be biased high, then they will 
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Table 2. For each luminosity bin, the number of lens galaxies, 
(z), and cr[z) (a characteristic width, though the distribution is 
not Gaussian), the mean weighted redshift z^ff, and the mean 
weighted luminosity L^ff relative to L« . The mean weighted red- 
shifts differ for the different source samples, and the results shown 
are the average of those values. 



Sample, Mr range Ngai (z) a-{z) z^g L^g/L, 



Blue lenses 



L3, 


-20, 


-19) 


215 808 


0.12 


0.04 


0.12 


0.5 


L4, 


[-21, 


-20) 


310 161 


0.19 


0.05 


0.17 


1.1 


L5, 


-22, 


-21) 


145 136 


0.27 


0.09 


0.23 


2.5 
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overestimate the value of r, and lead to an overestimate of 
/; if they are biased low, they will underestimate the value 
of r and consequently of /. Were it not for this problem, 
we would make a larger lens sample by going to fainter 
magnitudes, but the bias and scatter in the photometric 
redshifts at fainter magnitudes makes this impractical. As 
shown in M05 using data from DEEP2, for r < 19, the pho- 
tometric redshifts are not noticably biased, and have scatter 
Az ^ 0.04. Since Da{z) is not linear in redshift and / is not 
linear in r, in principle even with no bias and a scatter of 
Az ^ 0.04, the errors in r may bias our results, since we use 
r = 0Da{z)(1 + z). However, calculations indicate that this 
potential bias is well below the statistical error. 

The area covered by the full catalog is roughly 6200 
square degrees. Shape measurements were obtained for 
roughly 96 per cent of lenses passing the magnitude cut in 
this region, where many of the failures were due to saturated 
centres or other problems for very bright galaxies. 

The lenses were split into luminosity bins one model 
magnitude wide; the notation for the bins here is the 
same as in M05, but we only use the four brightest bins 
from that paper, for which the signal was detected with 
high signal-to-noise. The luminosities were computed us- 
ing the photometric redshift to get the dist ance modulus, 
and k- corrections are from kcorrect v1_11 iBlanton et alJ 
liooli). As in M05, a lum i nosity evolution correction con- 
sistent with lBlanton etlil] (|2003lJ) was applied, shifting all 
r-band magnitudes by -1-1. 6(z — 0.1). Table |5| includes in- 
formation about the luminosity bins, including the numbers 
of lenses and parameters of the redshift and magnitude dis- 
tribution. The effective redshifts and luminosities are com- 
puted using the same weights as are used in the computa- 
tion of the lensing signal. Information is shown for red and 
blue samples separately, where colour separation will be de- 
scribed shortly. 

In addition to the real lens catalog, our analysis also 
requires catalogs of random lens positions. These were cre- 
ated by distributing the random lenses uniformly across the 
survey area, since the photometric survey at r < 19 is fairly 
uniform and the density of sources at bright magnitudes 
does not vary significantly with seeing (at fainter magni- 
tudes, where star/galaxy separation is more difficult, the 
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Figure 3. The distribution of f rac.deV values for red and blue 

samples of galaxies, respectively 

density of sources does not show the expected increase with 
magnitude, and shows significant variation with seeing). 

It is also necessary to split the lens sample approxi- 
mately into elliptical and spiral samples. For this purpose, 
we used the obse rved u — r (model) colour, with a divi- 
sion at 2.22 as in lStrateva et all (1200 ll) . As a check of this 
colour separator, we used an independent method of classify- 
ing morphology, the frac_deV parameter output by Photo 
which, for each band, is determined by fitting the profile 
to frac_deV times the best-fit deVaucouleurs profile plus 
1— frac_deV times the best-fit exponential profile, and trun- 
cating the allowed range of values at and 1. This procedure 
is done in each band; we use the r band result due to its high 
signal to noise. Thus, frac_deV is determined only via the 
profile shape in one band, not from the colour. Fig.|3]shows 
the distribution of f rac_deV values for the red versus for the 
blue galaxies; as shown, the distribution for blue galaxies 
is strongly peaked at zero and for red galaxies at one, as 
would be expected if red colours and deVaucouleurs profiles 
are correlated, and if blue colours and exponential profiles 
are correlated. Hence we conclude that our "blue" and "red" 
samples include relatively pure samples (at the ~ 90 per cent 
level) of spiral and elliptical galaxies respectively. 

3.2 Source sample 

The source catalog is the same as in M05. In brief, shape 
measurements were performed on all SDSS imaging data col- 
lected through June 2004, using the re-Gaussianization tech- 
nique described in £irata & Scliak ( 20Qi), yielding roughly 
3 X 10^ source galaxies. Star-galaxy separation was achieved 
using Photo output OBJC_TYPE, and galaxies were re- 
quired to meet the following criteria: no flags from Photo 
indicating problems with object measurement (e.g., a satu- 
rated centre) can be set, r band extinction less than 0.2, r 
band extinction-corrected model magnitude r < 21.8, res- 
olution factor R2 > 1/3 in both bands (r and i) used for 
the shape measurement to ensure high-quality shape mea- 
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surements and to eliminate stars that were mis-classified 
by OBJC_TYPE as galaxies. For the purpose of the cal- 
culation of the shear signal, there are 3 disjoint sets of 
sources: "bright" (r < 21) sources, "faint" (21 < r < 21.8) 
sources, and high-redshift (0.4 < z < 0.65) Luminous Red 
Galaxies, or LRGs, which are explicitly excluded from the 
other two samples. The shear calibration was found in M05 
to be accurate to within [—7, -1-12] per cent {2a level) for 
r < 21 sources, [— 10, -fl8] per cent for r > 21 sources, 
and [—8, 19] per cent for high-redshift LRGs. Redshift dis- 
tributions were determ ined for faint (r > 21) sources using 
DEE P2 spectroscopy ("Dav is et. al I2003L iMadgwick et all 
I2OO3I. [Davis ct al. 2004, Co il et al.ll2004ll . whereas redshift 
estimates for individual galaxies were obtained for the other 
two source samples. For the bright (r < 21) sources, photo - 
metric redshifts from kphotoz v3_2 llBlanton et al.ll2003al) 
with error distributions determined using DEEP2 spec- 
troscopy in M05 were use d; LRG photometric r e dshift s with 
error distributions from IPadmanabhan et. al.l ll2005l) were 
used for the LRG sample. 



3.3 Systematics 

In M05, a large number of systematic errors in the weak lens- 
ing signal were analyzed to assess their significance in the 
SDSS data. For this work, many of the calibration uncertain- 
ties considered there (e.g., shear calibration bias, redshift 
distribution systematics, and stellar contamination) are ir- 
relevant due to our interest in the ratio of /AE to AS, both 
of which have the same calibration biases. However, there 
are other, more subtle systematics that are important here, 
that we describe in the subsections that follow. 



3.3.1 Systematic shear 

One possible contaminant of the / measurement is system- 
atic shear. As discussed in M05, any average smearing of 
the images along the scan direction will cause a slight align- 
ment of the lens and source ellipticities. This systematic thus 
raises (lowers) the shear signal on large scales by an additive 
factor if the PSF correction scheme over- (under-)estimates 
the correction to the shears due to this smearing. A constant 
systematic shear only changes AE on large scales, because 
for an isotropic distribution of sources around a lens, all 
contributions from the systematic shear to the tangential 
ellipticities cancel out; consequently, it is only important 
for lenses with an anisotropic distribution of sources. For- 
tunately, as described there, we can correct for the effects 
of systematic shear in the lensing signal using random lens 
catalogs to compute the signal due to systematic shear, so 
that we can subtract it off from the real signal. 

Unfortunately, for the /AE measurement, this proce- 
dure is insufficient. The alignment of lens and source elliptic- 
ities due to systematic shear will lower /AE (since 'jt.sys < 
where cos {2A6) > 0, and vice versa) thus lowering our esti- 
mate of /. This signal cannot be calculated analytically and 
removed without a very well-understood model for the size of 
the systematic shear and its variation on the sky (which we 
lack). However, there is one way to measure it: the system- 
atic shear contributes equally to /AE and to /45AS (since 
rotating the lens and source ellipticities by the same angle 



does not change the fact that they are aligned). Since the 
spurious shear contributions to /45AE and /AE are equal, 
we can subtract /45AE from /AE to measure the uncon- 
taminated value of (/ — /45)AEi3o, for which we presented 
model predictions in t|2.3l As theoretical justification for the 
subtraction to eliminate systematics, we note that the shear 
correlation functions are related via 



^++(e)±c.4^) 



1 

2^ 



ede [PE{e) ± PbW] 



(22) 

Here, ^++ and £,xx, the correlation functions of tangential 
and 45-degree rotated shear, are related via multiplicative 
factors to the measured /AE and /45AE45 (respectively) 
resulting from systematic shear, Pe is the power spectrum 
of the -E-mode contribution to systematic shear, and Pb is 
the power spectrum of the B-mode systematic shear. On 
small scales {£6 <^ 1, so Jo ^ 1 and J4 — > 0), ^++ = (,xx, 
i.e. the contribution of a constant systematic shear to /AE 
and /45AE is identical. As we will see, the signal from the 
systematic shear is indeed consistent with being associated 
with scales larger than those we use for the halo elliptic- 
ity measurement, justifying our scheme for eliminating this 
systematic. 



3.3.2 Anisotropic source number density 

In this paper, we are primarily concerned with transverse 
separations from 20 ft~^kpc (the minimum scale at which 
signal is measured) to several hundred fe~'^kpc, around the 
virial radius. For our measurement, we use a boost factor 
B{r) (to account for dilution of the signal due to inclu- 
sion of non-lensed galaxies in the source sample) that is 
azimuthally averaged, and therefore does not include vari- 
ations of the observed galaxy number density with AO. We 
consider here four effects that may cause azimuthal variation 
of the number density of galaxies around lenses, and their 
implications for our measurement of halo ellipticity using an 
azimuthally-averaged B{r): anisotropic magnification bias, 
sky subtraction errors, density-shape intrinsic alignments, 
and the inclusion of foregrounds in the source sample. 

The first such concern that can cause azimuthal varia- 
tion of the source galaxy density is magnification bias, the 
effect of which is determined for each source sample as de- 
scribed in section 4.7.3 of M05. The magnification bias on 
average affects the amplitude of both AE and /AE by an 
overall factor (since we boost by 1 + determined us- 
ing random catalogs, and therefore we are assuming that 
some of the real sources visible due solely to magnification 
bias are not lensed, so we overestimate the average signal). 
However, since magnification bias causes 5N/N oc k, and 
K varies azimuthally in our model, we must ask whether 
our use of azimuthally-averaged B{r) biases the measure- 
ment of fh. If the ellipticities of light and dark matter are 
aligned, then the magnification bias anisotropy will increase 
the source density along the lens major axis more than along 
the minor axis. This effect can be modeled as an isotropic 
background of sources that would be seen in the absence of 
lensing, plus an isotropic distribution of physically associ- 
ated sources, which together give us a boost factor Biso{r). 
Magnification bias leads to some additional boost factor, 
Bmaa{r)[l + 2 fmagCg cos {2A6)]. Howcver, in reality these 
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Figure 4. fmagir) for the NFW and TIS models as a function 
of of transverse separation, in units of scale radii. 



extra sources along the major axis are truly lensed, so there 
is no need for our boost factor to have azimuthal variation to 
account for the azimuthal variation in source number den- 
sity. Therefore, in the presence of anisotropic magnification 
bias, our use of the angle- averaged B{r) = Biao{r)+Bmag{r) 
will overestimate the lensing signal, biasing both AE and 
/AE, but not fh. 

While anisotropy of the source distribution around 
lenses due to magnification bias does not bias the measure- 
ment of / using an azimuthally-averaged boost factor, it 
does provide an interesting test of our results if we can de- 
tect it, because the values of / and fmag are related to fh in 
a consistent way by the halo density profile p(r). As shown 
in Eq.|51 for a purely power-law profile, fmag = fhCi/4:. Fig. 3] 
shows fmag as a function of scale for NFW and TIS profiles. 
When studying our results in i|K|we will return to the issue 
of detecting anisotropy of magnification bias. 

The second effect we discuss that may cause number 
density anisotropy is small-scale errors due to Photo, the 
sky subtraction systematic described in section 6.3.7 of M05. 
In short, a problem with determination of the sky level near 
bright (r <~ 17.5) lenses leads to variation of galaxy den- 
sity within about 1' due to errors in the sky flux scatter- 
ing galaxies into or out of our flux and apparent size cuts. 
Furthermore, model fits to determine fluxes and shapes are 
influenced in a difficult to quantify way. Consequently, we 
cannot assume an error in the amplitude of AE due to this 
effect that would cancel out of the ratio (/AE)/AE, we must 
ask whether it will induce alignments of densities on small 
scales that would change /AE but not AE. Furthermore, 
we must consider that this effect depends on the angle from 
the lens major axis, since at fixed transverse separations, 
the light on the major axis is brighter than on the minor 
axis. Fortunately, the analysis is simplified by the fact that 
in our r < 19 lens sample used for this work, only about 17 
per cent are at r < 17.5 where this effect is a problem. 

We may naively expect that if the sky systematic acts 
only affects number counts in our catalog, then they will be 



changed more on the major axis than on the minor axis. 
However, this systematic affects both lensed and physically- 
associated galaxies in the source sample, so dilution due 
to physically-associated sources should not have azimuthal 
dependence. Consequently, our use of azimuthally-averaged 
B{r) is indeed appropriate. 

The third effect we consider is the effect of intrinsic 
alignments. There are two intrinsic alignment effects that 
must be taken into account; the first is a correlation between 
the shapes of physically-associated galaxies, or the "shape- 
shape" correlation; the second effect is a correlation between 
the shapes of galaxies with the local density field, or the 
"density-shape" correlation. In this section we focus on the 
"density-shape" correlation, since it modulates the observed 
galaxy density around the lenses, and defer discussion of the 
"shape-shape" correlation for the next section. 

While we attempt to eliminate physically-associated 
sources from our calculation using photometric redshifts for 
the r < 21 sample, the fact that our boost factors are greater 
than one indicates that we are not fully successful. If the 
distribution of physically-associated sources around lenses 
is not isotropic, but is instead aligned with the lens ma- 
jor axis, then the number density of "sources" on the sky 
may also appear to have such an alignment due to con- 
tamination. Such an effect is caused by an alignment be- 
tween the lens galaxy shape and the local density field. This 
anisotropy of the satellite distribution has actually been de- 
tected before, albeit wi th confiicting am plitude and sign. 
The initial detection by lHolmber3 il969l) , which gave rise 
to the name "the Holmberg effect," was of an alignment of 
satellite galaxies along the minor axis of disk galaxies to 
transverse separations of 50 k pc. Since then, several stud- 
ies have shown a similar effect: IZaritskv et alJ lll997l) found 
alignment of satellites along the min or axis of spiral galaxies 
for separations of 300-500 kpc, and lSales fc Lambad ll2004ll 
found the same using a mixed sample of satellites in the 
2-degree Field survey (2dF) out to 500 kpc. Other stud - 
ies have shown the opposite; for example. iBrainerdI i2005l) . 
using a sample of satellites in the SDSS to projected sepa- 
rations of 400 fe~^kpc, fo und an alignment along the major 
axis; Man delbaum et al.t ( 2006) found a larger-scale align- 
ment (0.3-60 /i~^Mpc) along the major axis for L > L» 
spectroscopic galaxies in the SDSS. A nd finally, o ther works 
(Hawlcv & Peebles 1975 and MacGilli vrav et al.l 1982) have 
found isotropic satellite distributions. These results were 
found using small primary-satellite systems; results with 
clusters are more definitive, indicating a tendency for the 
bright cluster galaxy (BCG) to be elongated along the clus- 
ter major axis, and hence for it to have satellites concen- 
trated along /S.9 = 0, tt and a defi cit at 7r/2 and 37r/2, a man- 
ifestation of the Binggeli effec t ||Binggel| 198j: iFuUer et aD 
119991: IWest fc Blakesled I2OOOI: iKim et all 12003) . An align- 
ment of physically-associated galaxies with the lens major 
or minor axis poses a problem for our measurement of /: if 
there is a larger (smaller) number of physically associated 
"sources" along the major axis, our azimuthally-averaged 
boost factor underestimates (overestimates) the signal on 
the major axis and vice versa for the minor axis, thereby 
underestimating (overestimating) /. We attempt to mea- 
sure this effect using a test to be described at the end of 
this section. 

The final effect we consider is the inclusion of galaxies 
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in the source sample that are actually between us and the 
"lens," and that lens the "lens" (foreground contamination) . 
Because this effect causes the "lenses" to have tangential el- 
lipticity relative to the "sources," it may increase number 
counts of "sources" along the "lens" minor axis, and de- 
creased number counts along the major axis. Because these 
foregrounds are not lensed, our boost factor needs to have 
azimuthal dependence to account for them; since it does not, 
we may overestimate /. (They also lead to overall suppres- 
sion of the signal which is unimportant for our / measure- 
ment.) Due to our cuts, we do not expect there to be large 
numbers of foregrounds in our sample, so this effect should 
not be a major problem, but still test our assumption using 
the tests described below. 

To determine which, if any, of these four possible causes 
of azimuthally-dependent variations in source density actu- 
ally occurs, we use the following test. For each lens sample 
and the sources with photometric redshifts, we divided the 
lens-source pairs into three samples: the "foreground" sam- 
ple with Za < zi ~ e (e = 0.05, 0.1), the "physically associ- 
ated" sample with zi — e < Zs < zi + e, and the "source" 
sample with Zs > zi + e. Due to photometric redshift errors, 
these divisions are imperfect. For each of these three samples 
and lens luminosity /colour samples, we compute the angle- 
averaged boost-factor B{r) = Nls/Nrs, where Nls is the 
number of pairs of real lenses and real sources, and Nrs is 
the number of pairs of random lenses and real sources (with 
weight equal to 1 for each galaxy) , and we also compute the 
sum over pairs ^Ae(r) — X^ls '"^^ {2A6) /Nrs ■ 

If we assume that the galaxy density has some isotropic 
distribution Ni3o{r) plus a varying component A'^Ae('')[l + 
'^fsysSg cos (2A^)] duc to One of the aforementioned effects, 
then we will have 

B{r) = [iV..o(r) + iVAfl(r)]/iVfls(r) (23) 
Ueir) = 2NAe{r)fsysel^JNRs{r) 

(the factor of two comes from our use of the full elliptic- 
ity Bg, the sum of the squares of the two components, for 
which (eg) = 2erma)- Comparing ^A9{r) with B{r) allows 
us to determine the effective feff{r) of the distribution of 
sources, where feff{r) = fsysNAe{r) /{NAe{r) + Niso{r)) (if 
the anisotropic component of the number density is much 
smaller than the total angle-averaged number density, the 
effect of the anisotropy will be correspondingly smaller). 
For the two sources of anisotropy NAe{r) for which our use 
of azimuthally-averaged boost factor is incorrect (intrinsic 
alignments and foregrounds), we will have underestimated 
/ and therefore / — /45 by the feff corresponding to that 
effect, to lowest order in e. This underestimate can be un- 
derstood considering that AE cx (1 -I- 2/eg cos (2AS)), but 
if the boost factor has azimuthal variation, then we really 
measure shear oc (1 + 2f-yeg cos (2A6')) such that when mul- 
tiplied by the boost (x (1 -I- 2/e//eg cos (2A6')), we get the 
true AE (/e// -I- /7 ~ f)- When we neglect the /e// of the 
boost, we therefore underestimate / by this additive factor. 
However, when we measure /45AE45 oc /45eg cos (2A^), to 
lowest order in eg there is no correction to account for our 
use of the azimuthally-averaged boost. 

We expect a variation in number density due to mag- 
nification bias will appear predominantly on small scales in 
the "source" sample (most strongly for bright lenses); due to 
intrinsic alignments, mainly in the "physically associated" 



sample but with uncertain scale dependence; due to Photo 
effects in all three source samples equally, but only on small 
scales; and due to lensing by foregrounds predominantly in 
the "foreground" sample. Results of these tests will be pre- 
sented in H5.ll 

3.3.3 Shape-shape intrinsic alignments 

As mentioned in >I3.3.2I physically-assocated lens-source 
pairs may potentially contaminate our measurement due to 
any shape-shape intrinsic alignments. If the degree of galaxy 
alignment is independent of A^, then the shape-shape align- 
ment manifests in the same way as a constant systematic 
shear, and will be eliminated by the /AE — /45AE45 sub- 
traction. However, any azimuthal structure in the shape- 
shape intrinsic alignments would not be removed via this 
subtraction, and will contaminate the measurement. To test 
for this problem, we calculate /7 — /is 745 for sources in the 
"physically associated" sample (we measure this instead of 
/AE — /45 AE45 because E^^ is zero for Zs < zi). Results of 
this test are presented in section 

5.5.^ Y'nO'TO-related shape measurement biases 

In principle, the sky brightness systematic may bias the 
shapes of sources in a way that would affect the measure- 
ment of /AE — /45AE45. This might occur if the sources 
align with the gradient of the sky brightness near the lenses. 

To test for this contamination, we measure /7 — /45745 
using sources in the "foreground" sample (with Zs < zi ~ 
0.1). This quantity should not have significant contami- 
nation from actual lensing or from azimuthal structure in 
intrinsic alignments, and hence should be dominated by 
software-related systematics. Results of this test are pre- 
sented in section 1^31 



4 METHOD OF ANALYSIS 

The first step for the halo ellipticity calculation is to com- 
pute AE, /AE, and /45AE as in Eqs. |S|and refE:sums45 
and the accompanying text. As described in tl3.3.1l /45AE 
will allow us to identify systematic shear contributions to 
/AE. Covariance matrices for these quantities are computed 
via bootstrap resampling with 150 subregions, similar to the 
method used in M05. Once we compute these quantities as 
a function of comoving pair separation r, there are several 
ways to proceed with the analysis. 

The first method used is a non-parametric determina- 
tion of / — /45 averaged over radius. We start with the signal 
AE for the desired range of radii, and /AE — /45AE with 
errors determined from those quantities added in quadra- 
ture (they are uncorrelated). We then consider the prob- 
lem of two Gaussian variables y and x that are related via 
m — y/x (here y represents /AE — /45AE, x represents 
AE, and m is our desired quantity, / — /45). In our case, 
we have multiple estimators yi and Xi (the values at each 
radial bin) and would like to combine them while takiiig 
into account non-Gaussi anity in the ratio y/x |Bliss 1935a| . 
lBlissll935bl . lFiellerll954^ . For each measurement i, the quan- 
tity yi — mxi is a random Gaussian variable drawn from an 
A''(0, (Jy. -l-m^a^.) = N(Q,w~^) distribution. Consequently, 



© 0000 RAS, MNRAS 000, 000-000 



10 Mandelbaum et. al. 



the following summation over all measurements is also a 
random Gaussian variable: 



X] Wi(yi - mxj 



N 



1 



^ , (24) 

where the distribution is taken at fixed m, and the weights 
Wi depend on m. We can then determine frequentist confi- 
dence intervals^ at the Za level by writing 



J2 Wijiji — mxj) 



< 



(25) 



Since this expression depends on m, in order to calculate the 
confidence intervals we make a grid in m and interpolate 
to find the value of m for the desired confidence interval, 
e.g. m{Z = 0) corresponds to the average value, m(Z = 
±1) give the 68 per cent confidence limits, and so on. The 
value of this approach is that it gives a radius- averaged value 
of / — /45 without requiring a model for the shape of the 
profile or of / — f^s, that it easily allows us to combine 
measurements of /— /45 with different source samples for the 
same lenses, and that it takes into account non-Gaussianity 
of the error distributions of the ratios yi/xi. Once we have 
measured {/ — /45)r nonparametrically, we relate the result 
to the predictions for different profiles to extract the value 
of interest, f^. 

The second method used here is to fit the signal to a 
power-law profile using the predicted relationships in Eg. 1141 
That is, we write 



AE(r) = Ar~ 
/AE(r) - /45AE(r) = A,h 



(26) 



(a -2)2 



4a 



and fit jointly for A, a, and fh- Unlike the previous method, 
it does not take into account non-Gaussianity in the error 
distribution of fh, which we will see is important (and we 
cannot use this method to measure fh with the NFW or TIS 
density profiles). 

Finally, we repeat this analysis on large scales (800 
/i~^kpc to 2 /i~^Mpc) with all samples to check for align- 
ment of the light distributions of spirals and ellipticals with 
local LSS. 



5 RESULTS 

5.1 Number density systematics 

Here we present results related to the azimuthally- 
dependent boost factor tests described in tl3.3l As de- 
scribed there, for "foreground," "physically associated," and 
"source" samples (determined using photometric redshifts) , 
we compute the effective feff{r) with each of the lens sam- 
ples; this computation allows us to determine whether mag- 
nification bias, Photo effects, intrinsic alignments, and lens- 
ing by intervening foregrounds may be significant. 

For blue lenses of all luminosities, and for L3 and L4 red 
lenses, fsff{r) was found to be consistent with zero on all 
scales, as will be discussed further. For red L5 and L6 lenses, 
this was not the case; a plot of feffir) and B{r) for these 
two lens samples is shown in Fig. O with the three source 
samples. 

^ We are not constructing a Bayesian posterior region P{m). 
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Figure 5. The measured feffix) and B{r) for red lenses, L5 and 
L6, with the three different sets of sources: foregrounds, physically 
associated, and sources. Note the different scales on the vertical 
axes for L5 and L6. 



There are several important features in this plot. First, 
we consider the bottom two plots, B{r) for L5 and L6. The 
fraction of physically associated sources in a given sample is 
equal to {B — 1)/B. So, we see that our photometric redshift 
cuts are reasonably efficient at isolating a sample of phys- 
ically associated sources; at 100 /i~^kpc, in L6, nearly 2/3 
of the "physically associated" sample is actually physically 
associated, compared to about 15 per cent of the "source" 
sample and 35 per cent of the "foreground" sample. The 
fact that B{r) for the physically associated sample actually 
increases with radius for 20 < r < 30 /i~^kpc in L5 and 
20 < r < 40 /i~^kpc in L6 indicates that these ranges of 
radii are affected either by the sky subtraction problem, or 
by the actual light from the lens galaxy; in any case, due 
to the possibility of photometric errors in these regions, we 
exclude them from the analysis in the rest of the paper. 

Next, we consider the top plots in Fig.|3 feffir) for L5 
and L6 red lenses. For L6, we can immediately see that for 
100 < r < 2000 h~^kpc, feff{r) is positive for the "physi- 
cally associated" sample in a way that is highly statistically 
significant, but not for the other samples. This fact points 
to density-shape intrinsic alignments as the cause for this 
finding. Apparently, there is a tendency for satellite galax- 
ies in these ranges of radii to cluster along the major axis 
of the "lens" galaxy. The results shown here include every- 
thing from small groups to large clusters, so what we see 
is an average of this effect over systems of all sizes. An- 
other work (in preparation) will further explore the lumi- 
nosity and environment-dependence of this effect in an at- 
tempt to understand the cause. This plot shows feff, the 
actual value of fint associated with intrinsic alignments (as- 
suming that the number density on the sky about lenses is 
isotropic with the exception of this effect) is determined via 
feffir) = (B(r) - l)/,„t(r)/B(r). Fig.ilshows /,„i(r) for 
this sample. 
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Figure 6. The measured finti^), including non-Gaussian er- 
rorbars, for L6 red lenses derived from the results for physically 
associated sources 



We could also attempt to derive similar results from the 
foreground or source sample (at least for r ~ 100 /i~^kpc, 
where there is a slight signal) for the sake of comparison. 
Since B{r) is much lower, feffir) due to intrinsic align- 
ments is much lower for these samples. When the results 
are compared for these radii, the resulting fint is consistent 
with that shown in Fig.|S|at the la level (we may not expect 
them to a priori be identical, since the samples are different: 
the physically associated sample includes those with correct 
photometric redshifts, whereas the foreground and source 
sample includes physically associated sources with incorrect 
photometric redshifts, and since they lie in separate regions 
of colour space, they may have different properties with re- 
spect to intrinsic alignments). 

For L6, it is clear that on scales larger than 100 /!,~^kpc, 
the anisotropy of the number density distribution around 
lenses is due to density-shape intrinsic alignments. Now, we 
consider scales r < 40 /i~^kpc; we already know that the 
sky subtraction effect or the light of these lenses is caus- 
ing a loss of number density (as in the B{r) plot), and the 
rapid decrease in feff{r) with decreasing r on these scales 
indicates that the suppression of number density is stronger 
along the major axis of the lens light distribution, as may 
be expected. This explanation is supported by the fact that 
on these scales feff{r) has roughly the same pattern for 
all three source samples, and therefore the cause cannot be 
a physical effect on only foregrounds, physically associated 
galaxies, or lensed sources. This finding supports our deci- 
sion to not use these scales when attempting to measure 
halo ellipticity. 

Finally, we consider 40 < r < 100 h~^kpc. For this 
range of radii, ,feff{r) changes with r in similar ways for the 
foreground and physically associated samples; this suggests 
that it is due to physically associated sources in both, with 
the amplitude difference being due to the different fraction 
of physically associated sources in the two samples. However, 
as r decreases in this range, feff{r) increases much faster for 



the lensed sources; this fact suggests that the increase is due 
to anisotropy of the magnification bias causing an increase 
in source number density preferentially along the lens major 
axis. Unfortunately, the exact /e// for this effect is difficult 
to extract from the plot due to the unknown radial depen- 
dence of the feff due to intrinsic alignments. We could also 
attempt to predict it, using the assumed value of fmag{r) 
for a given profile, with the appropriate 5N/N, to get 

/ ^ , .SN/N 

feff,mag{r) = fmag{r)^^ (27) 

(The division by B{r) is because we need the ratio of the 
additional lensed sources to the total number of pairs, not 
to the total number of lensed pairs which is given by SN/N.) 
The difficulty is that it is hard to untangle the effects of fint 
and fmag in a model-independent way, which means that (1) 
we cannot easily use the measured fmag as a check on the 
results for /, and (2) we cannot easily correct the measured 
/ for the effects of fi„t (since it, but not fmag, contaminates 
the measurement of fh)- 

We can consider several approaches to this problem, but 
the simplest is to assume that the intrinsic alignments in all 
three source samples are similar (i.e. there is no colour or 
magnitude dependence), and therefore to get feff.intij) for 
the lensed sources, we can just use the implied fint{r) from 
the "physically-associated" sample in Fig.|S|with B{r) to get 
feff,int{r) for any source sample. Comparing this prediction 
against the measured feff{r) for a given source sample al- 
lows us to untangle the effects of fint and fmag- 

For L5, the situation is similar, though (1) the small- 
scale suppression of number counts only extends to 30 
/i~^kpc, not 40 /i~^kpc, and (2) the detection of ellipticity 
of the satellite galaxy distribution (density-shape alignment) 
implied by the nonzero f^f / (r) for the physically associated 
sample is not as large as for L6, though still statistically sig- 
nificant for r < 400 /i~^kpc. There is a hint of magnification 
bias anisotropy on scales up to 70 /i^^kpc for this luminosity 
bin. 

The results shown in Fig. [pimply that we do not observe 
significant lensing by intervening foregrounds that may have 
biased our fh measurement. This may be because the SDSS 
data are rather shallow. 

We now ask how these findings will affect the results 
with fainter sources as well. There are three effects that we 
have observed: small-scale effects due to sky subtraction or 
the lens light distribution causing a loss of sources prefer- 
entially along the major axis; effects on all scales due to 
contamination by physically associated sources that tend to 
align along the lens major axis; and anisotropic magnifica- 
tion bias for r < 100 /i~^kpc (L6 red lenses), which increases 
source number density preferentially along the lens major 
axis. (Since foreground contamination is not a problem for 
the r < 21 sources, at the lowest redshifts, they are cer- 
tainly not a problem for the fainter, higher redshift source 
samples.) We consider all three effects separately. 

First, we expect that the fainter sources will be affected 
by the small-scale light problems similarly to the r < 21 
sources used in this section. Consequently, for these sources 
we also avoid r < 30 /i~^kpc (L5 red lenses) and r < 40 
/i~^kpc (L6 red lenses), but use all scales above 20 /i^^kpc 
for other lens samples. 

Second, the feff{r) due to physically associated sources 
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is the product of two factors, fint{r) and {B{r) — 1)/B{r). 
For the fainter sources, B{r) is about 30 per cent higher for 
the fainter sources (due to the lack of photometric redshift 
information) than for the brighter ones for r > 60 /i~^kpc; 
below that separation, the difference is not as large. For 
small boosts as in all luminosity bins for r > 100 h'^kpc and 
even for smaller scales for L3 and L4, this translates roughly 
to a doubling of the {B{r) — 1)/B{r) factor used to determine 
the feff. However, it is not clear that the fint{r) for these 
sources should be the same as for the brighter sources. To 
determine whether this is a safe assumption, we repeated 
the tests shown above with r < 20.3 and 20.3 < r < 21 
sources separately to check for evolution of fi„t with appar- 
ent magnitude. In both L5 and L6, while feffir) and B{r) 
were different for the r < 20.3 and r > 20.3 sources, the im- 
plied firit(r) values were consistent with each other at the 
la level for all ranges of radii and all luminosity samples 
that we use for this analysis. We thus use the same fint{r) 
for r > 21 sources as for r < 21 sources; if this is not the 
case, we will see that the correction we apply is quite small 
compared to our statistical errorbars. For LRG sources, the 
lens and source samples are essentially uncorrelated, so con- 
tamination by physically associated sources is not a concern. 

Finally, the magnification bias effect depends on two 
things: the magnitude distribution of the source sample, and 
the R2 distribution. In M05, we derived that the predicted 
5N/N due to magnification bias is 1.9k, 0.7k, and 3.1k for 
r < 21, r > 21, and LRG sources in our catalog, respectively. 
K oc E~^, which also varies with source sample. 

When we actually use the data, we must use the 
weighted (i.e., with weights used for lensing) B{r) and 
feff{r) when doing the corrections. The weighted B{r) for 
the data can be significantly lower than the unweighted, 
since the E^^ weighting weights higher-redshift sources more 
heavily, and those sources are less likely to be physically as- 
sociated with the lenses. 

We want to use the weighted B{r) and the angular 
correlations to determine fef f,int{r), the portion of the 
anisotropy in the source distribution that is due to contam- 
ination by anisotropically-distributed physically-associated 
pairs, excluding the portion of feff that is due to magnifi- 
cation bias. We then add feff,int to the measured / — /45 
to find the non-contaminated value. The correction we cal- 
culate is well below the la errors, so too much importance 
should not be placed on details of the correction scheme, 
which we outline briefiy here. We assume that fint is inde- 
pendent of source sample, and fmag for each sample must 
be related by consistency relations due to its dependence 
on the invariant E. We fit the fint value to an empirical 
model (linear in log(r), which has no theoretical meaning) 
and extrapolate the values of fint to lower radii to see what 
feff,int they predict, assuming that the residual f^jj is due 
to magnification bias. 

The resulting feff,int{r), the amount by which we un- 
derestimate the measured / — /45 due to our use of the 
azimuthally-averaged B{r), is then averaged over the range 
of radii used for the measurement. We find that f(r) — f4,5{r) 
is underestimated by 0.01 for r < 21 and r > 21 sources with 
L5 red lenses, and by 0.02 for r < 21 and 0.03 for r > 21 
sources with L6 red lenses. Intrinsic alignments are not a 
significant contaminant for the LRG sources because they 
are at higher redshifts than the lenses. All of the results in 



AE, fAE, and f„AE for L4 red lenses 
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Figure 7. Plot of measured AS (top), /AS, and /45AS (bot- 
tom) for red L4 lenses with r < 21 sources. 

the following sections are automatically corrected for this 
effect. 

5.2 Systematic shear 

Here we present results for AE, /AE, and /45AE for the 4 
luminosity bins and 2 colour samples. As an example of our 
findings, Fig.|7|shows these three quantities for L4 red lenses 
with r < 21 sources. As shown, /45AE is clearly negative 
at small transverse separations, and becomes closer to zero 
with increasing radius. 

We must then ask what could be causing this sig- 
nal. There are three possibilities to consider. First, we ask 
if cosmic shear could be causing this result (i.e. lensing 
by LSS at redshifts lower than the lens redshifts, caus- 
ing real alignment of lens and source ellipticities), so that 
the Csys that we are detecting is an actual effect rather 
than an artifact of some aspect of the shape measure- 
ment. We computed the expected cosmic shear signal for 
the combinations of lens and source samples with the high- 
est expected cosmic shear signal, assuming the best fit 6- 
parameter ACDM cos mologv from | Selj ak et alJ i2005lJ) . the 
transfer funct ion from [E isenstei n'fcHurin^Q^ . the growth 
function from Carroll et al. '("199^), and the non-linear map- 
ping from Peacock & Dodds (199^). The results in all cases 
indicate that cosmic shear is more than 10 times too small 
to cause the observed lens-source alignment. 

Another possibility to consider is that the spurious el- 
lipticity Esys has been imparted to the lenses and sources in 
the course of PSF correction. The PSF is determined in a 
given frame (10' by 14', or 1489 by 2048 pixels) using ±2 
frames (5 frames total) to determine the Karhunen-Loeve 
(KL) eigenfunctions, then using ±1 frame (3 total) to do a 
second-order polynomial fit to the coefficient s. With 15-25 
stars per frame used to determine the PSF iLupton et alJ 
l200ll) . the average separation between stars used to deter- 
mine the PSF is ~ 2.6'. There are a few ways that incorrect 
PSF determination could cause this effect. The first is sta- 
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tistical: stars from r ~ 16 to 19 are used to determine the 
PSF. In our source catalog, the error on the ellipticity for 
a galaxy at r ~ 19 is roughly 0.02. If even half of the stars 
used to determine the PSF are in the range 18 < r < 19, 
the statistical error on the PSF ellipticity is of the correct 
order to cause the effect we see. The spurious ellipticity Caya 
due to this effect will of course vary randomly from star to 
star, but since all stars are used to determine the PSF vari- 
ation across a frame, the statistical error due to a single star 
may affect a region larger than the separation between stars 
(since 9psf ~ 300 /i~^kpc at the mean effective redshift of 
the sample in Fig0. 

PSF systematics may also cause this effect; the PSF de- 
termination is accurate to within 2 per cent, so systematics 
in the PSF determination could give rise to some ellipticity 
of 0.01-0.02. These systematics include the possibility that 
the spatial change in the PSF is not really quadratic, or 
strange effects near the edges of the chip. 

Now that we know that /45AE and /AE are contam- 
inated in similar ways on small scales, our approach must 
be to use /AE — /45AE, and compare against model pre- 
dictions for this quantity on small scales where 6 <g Ocorr, 
of order a few hundred h~^kpc. 

As a further test of the systematic shear subtraction 
technique, we measured the correlation function of the sys- 
tematic shear from the data. To do so, we used data from 
the southern galactic regions, for which there were multiple 
(> 20) observations of the same area. For 30 pairs of runs in 
this region (selected randomly), we isolated the set of objects 
with shape measurements in both runs. For those objects, 
we then found pairs of galaxies as a function of angular sep- 
aration. In each run, we computed et and for each galaxy 
in the pair relative to the other, and constructed differences 
Aet and Ae^ for each pair, where the A signifies the differ- 
ence between the values of et and ej, in the two runs. Thus, 
Aet and Ae^c encode information about systematic differ- 
ences in the ellipticities of the same galaxy from run to run. 
For each pair, we then used the values of Aet and Ae^c to 
construct correlation functions 

CVT'W = (^et,iAet,2) (28) 

Cir'(e) = {Ae.,iAe,,2) (29) 

Cir>(e) = {Ae,,iAe..2> (30) 

where 1 and 2 denote galaxies in the pair. 

Fig. |H| shows the results of this computation, divided 
by two to approximate the result for the real data, because 
when comparing the ellipticities across runs, the system- 
atic shear in each of them presumably adds coherently. As 
shown, and ^i'^''"' are consistent with each other, and 

C+x"^ is consistent with zero, as expected for any parity- 
preserving systematic such as spurious ellipticity along the 
scan direction. Furthermore, the decline of the correlation 
functions with angular scale indicates that we are correct 
in our understanding that this is a small-scale systematic. 
For reference, we note that in L3, the typical radial range 
used for our calculations (20-280 h~^kpc) corresponds to 
0.0055 < 6* < 0.08 degrees, and in L6, 40-320 h'^pc corre- 
sponds to 0.0040 < 9 < 0.03 degrees. We see from the figure 
that it is possible that {Ae^Aex} is lower than (Ae+Ae+) 
by roughly 5 x 10~^, and have verified that the resulting 
changes in the computed values of fh are significantly less 
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Figure 8. Plot of measured correlation functions of the system- 
atic shear as described in the text. 

than half of the la errors for all models, where our over- 
estimate of (Ae^Aex) translates to an overestimate of fh- 



5.3 Shape-shape intrinsic alignments 

Here we present results of the calculation of f'y — /45745 for 
the "physically associated" source sample. We then compare 
the value of ,f"f — f 45^)4,5 to the observed AE to see what con- 
straints we can place on the contamination of fh from az- 
imuthal structure in shape-shape intrinsic alignments, fh'^ ■ 

The resulting f'y — 745745 is consistent with zero at the 
Icr level for all color and magnitude lens samples, on all 
scales. This result suggests that the contamination of fh 
must be consistent with zero as well; nonetheless, since the 
predicted values of / — /45 for some profiles are so small, we 
cannot formally place strong constraints on the contamina- 
tion. For example, for the red lenses, averaged over luminos- 
ity, the constraint on /^^^^ is 0.3 ± 0.6. For the individual 
luminosity bins, again for the red lenses, the resulting values 
of contamination are fi"^ = 0.0±1.8, 0.35±0.81, 0.57±1.0, 
and — 0.7 ± 1.7 for L3-L6 respect ively. 

We note that a recent work iMandelbaum et alj|200(?) . 
indicated no evidence of shape-shape intrinsic alignments on 
0.3-60 ft-^Mpc scales using a low-redshift {z ~ 0.1) sample 
of galaxies from the SDSS, despite the detection of a sta- 
tistically significant density-shape alignment. However, this 
measurement did not extend below the virial radius, and 
involved pairs of bright galaxies rather than one relatively 
bright and one significantly fainter, so we cannot rely solely 
on it. 



5.4 PHOTO-related biases 

In this section, we present results of the measurement of 
/7 ~ f 4,51 4,5 using the "foreground" sample selected using 
photometric redshifts. Any measurable /7 — 745745 should 
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Table 3. Central values and 68 per cent CL limits on the con- 
tamination of fh due to software biases. 



Luminosity sample f'f^^^ (blue) Z^*''*' (red) 



L3 3.7tg;5 1.3 ±3.5 

L4 -O-Stli -1.0 ±1.1 

L5 -0.6 ±2.4 -0.9 ±0.9 

L6 -0.4l^:? -0.4t°J 

All 0.0 ±1.8 -0.7 ±0.6 



thus be due to software-related biases, e.g. problems in de- 
termination of sky brightness or deblending errors. 

As for the previous section, we found that the measured 
/7 ~ /45745 for this sample was completely consistent with 
zero on all radii, for all lens samples. However, when we 
compared it against the real AE for the "source" sample, we 
found that we cannot place very strong constraints on the 
contamination of fh due to photometric errors, /^''"'.TableEl 
shows the constraints for each color and luminosity sample, 
computed using the SIS profile predictions for f — fis (results 
with an NFW or TIS model are qualitatively similar). 

5.5 Halo ellipticities 

Ifaving computed the signal, we proceed with a basic anal- 
ysis: a non-parametric determination of / — fis averaged 
over radial ranges for each lens sample and source sample 
combination, as in f0| There are several considerations for 
choosing radial ranges: first, we prefer scales small relative 
to 9corr, so that the /AE — /45 AE subtraction is justifiable; 
second, that it is not clear what signal to expect on scales for 
which neighboring galaxies in the same halo contribute sig- 
nificantly to the lensing signal, so we restrict to radial ranges 
in which the majority of the weak lensing signal is coming 
from the lens alone. For L3-L4, we start the measurement 
from the minimum pair separation used here, 20 /i~^kpc, 
and go as far as 245 h^^kpc (L3), and 300 /i~^kpc (L4, L5); 
for L5, we use 30-300 /i~^kpc, and for L6, we use 40-300 
h~^kpc because of the small-scale systematics noted in H5.ll 
The maximum separations are all smaller than 6corr, so the 
subtraction of /45AE accurately removes contamination due 
to systematic shear. Table |1| shows the average / — /45 for 
each lens and source sample combination. 

Several results in this table require comment. First, the 
/ — fi5 values using r < 21 sources appear more negative 
than those for the other source samples, so we must check 
for consistency of the samples. We do this by finding, for 
source samples i, 

(/)source = ^"^'^^2 (^1) 

then getting 

X' = J2 ( ^•^''^■^]"^' ) (32) 

which should follow a distribution with degrees of free- 
dom equal to one less than the number of samples being 
compared. We indeed find that the results with the difi'erent 
samples for each lens sample are statistically consistent. 



While the test is the most general one that is sen- 
sitive to all kinds of discrepancies, the fact that 5 out of 
the 8 values of fh for the r < 21 sources are lower than for 
the r > 21 sources suggests that we must also devise a test 
that is more sensitive to a systematic offset. For this test, we 
compute the eight values of Sf — — j(r>2i) ^^^^ their 

errors a{df), then find (Sf) via a weighted (by l/a^{Sf)) av- 
erage over 5f values, and the error on this quantity, a{{5f}). 
The result is that the weighted average (Sf) — —0.13, with 
an error <j{{Sf)) = 0.10. Consequently, the discrepancy is 
just over la, and therefore not statistically significant. 

As shown, the constraints for the blue galaxies are 
weaker than for the red galaxies. When we average over all 
source samples for a given luminosity bin for blue galaxies, 
we find results that are negative for L3-L5 but consistent 
with zero at the 1 — 1.5a level, and for L6, that are posi- 
tive but again consistent with zero at the la level. For red 
galaxies, the results are somewhat more complicated. The 
measured value of / — /45 appears to roughly increase with 
luminosity, which can be explained in several ways, as we 
will show. We do not show results averaged over lens lumi- 
nosities, because the results for / — /45 must first be related 
to predictions for / — /45 for a given profile to extract the 
measured fh for that model, and then averaged over luminos- 
ity. In other words, for a given profile, the resulting / — /45 
for these radial ranges are not the same for different lens 
luminosities even for the same fh (due to variation of Vs 
with luminosity), so an average of f — /45 over luminosity is 
meaningless. 

We now relate the non-parametric determination of 
/ — /45 in Table |1] to fh for various halo profiles. However, 
because our results for / — /45 have such large errors, we 
will not use an extremely detailed or precise model. The 
simplest model, a SIS, predicts f/fh = 0.25 and /45 — 
on all scales, so that for all the results above, we can write 
ff^ = - /4g)/0.25. We also consider the NFW and TIS 
models, but do not attempt to do detailed profile fitting, 
which is complicated by possible calibration errors due to 
use of photometric redshifts for foregrounds. What mat- 
ters is only an approximate value of rs for the NFW and 
TIS profiles. To get an approximate value of Va using the 
NFW profile, we use best-fit masses for these lens samples 
from th e hal o model fits as descr ibed in iMandelbaum et alJ 
^2005b^ and ISeliak et alJ ((20053). Since the cahbration for 
the photometric foreground samples is uncertain, we use the 
same luminosity and colour divisions with the spectroscopic 
sample lenses, and neglect evolution in the mass-luminosity 
relationship with redshift (though we do include passive lu- 
minosity evolution as noted in il3.1ll . For the given best-fit 
mass, we can compute the concentration parameter at the 
effective redshift of these lens samples using the relat ionship 
used in that work and in IMandelbaum et al.l ll2005lJ) . 

Furthermore, since for an NFW profile the mass and virial 
radius (where virial mass here is the mass within the radius 
at which the density is 180 times the mean density, or 54 
times the critical density with flm = 0.3) are related, the 
halo mass and redshift alone determine the scale radii for 
these luminosity bins. We then used the profile-dependent 
forms for / — f^s in Figs.0and|5|to redo the nonparametric 
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Table 4. Results for / — 745 averaged over radial ranges described in the text for each combination of lens and source samples. All 
errors shown are 68 per cent confidence intervals 



Lenses 


r < 21 sources 


r > 21 sources 

/-/45 


LRGs 


all 


L3, blue 
L4, blue 
L5, blue 
L6, blue 


-0.04 ±0.30 
-0 64+°"^ 
-0.01 ±0.43 

'-'■^'-0.80 


-0 42+0-43 

n 90 + 0.35 
^•^'-'-0.36 

0.7711° 


9 1+2.3 
^■^-3.0 

1.82tli 


— qR+0-25 
U.^>D_(j 26 

-0.27 ±0.28 

1 o+i^* 


L3, red 
L4, red 
L5, red 

L6, red 


-0.12 ±0.21 
-0.08 ±0.12 
-0.19 ±0.13 
0.04 ±0.22 


'^■^'^-0.21 
0.07 ±0.11 
0.12 ±0.12 
0.34 ±0.16 


-n 25+° 23 

'-'■^'-'-0.24 

0.04 ±0.15 
0.38 ±0.16 
0.41 ±0.24 


'^■'^■^-0.13 

0.01 ±0.07 
0.08 ± 0.08 
0.29 ±0.12 



determination of / — f4s described above dividing by (/ — 
fis)/ fh for our models in each bin, yielding a determination 
of fh for these two profiles that includes non-Gaussianity 
when averaging over these radial bins. For the TIS profile, 
we fit AE for the profiles to obtain scale radii only (since the 
amplitude is a nuisance parameter). Table |^ shows relevant 
model parameters for the SIS profiles, for the NFW model 
(Afi8o and r^), and for the TIS profile (rs), in addition to 
the resulting from the non-Gaussian averaging method. 
We note that the scale radii show the expected trend of 
increasing with luminosity , and that the best-fi t masses are 
consistent with those from lSeliak et all ll2005afl for each bin 
when we average over the results for each colour. 

First, we consider the results for blue lenses in Table |S| 
With the SIS profile, L3-L5 each have fh ~ If negative, 
and L6 has fh ~ la positive, yielding a net result of fh = 
— 1.1 ± 0.6, so almost 2a negative. We note that the results 
for each luminosity bin are statistically consistent at the la 
level, and that the result averaged over luminosity is not 
equal to that with 1/a^ weighting because it includes the 
non-Gaussianity of the error distributions. With the NFW 
profile, the measured fh values are larger in magnitude than 
with the SIS, as expected due to the lower predicted (/ — 
fis)/ fh for this model, and the averaged result is negative 
but consistent with zero; the same is true for the TIS. 

For red lenses, the assumption of an SIS profile gives a 
net value of fh = —0.06 ± 0.19, so consistent with no halo 
ellipticity. However, the NFW and TIS profiles, with lower 
if — fib)/ fh (even going slightly negative for the TIS), yield 
positive luminosity-averaged results that are quite similar to 
each other, 0.60 ± 0.38 (NFW) and 0.57 ± 0.41 (TIS), but 
again consistent with zero. We note that the errors for the 
NFW and TIS profiles are larger than for the SIS profile; this 
is likely due to the fact that the predicted ( / — /45 ) / fh values 
are so much lower and vary with radius, so since we divide 
by this small number on large scales, these scales do not 
dominate as much as they do for the SIS model, where the 
low measured values of fh tend to dominate the averaging 
process and lead to a small result with small errors. For red 
galaxies, the results in different luminosity bins appear to 
be statistically inconsistent; use of the test for sample 
consistency in Eq. [^yields = i for 3 degrees of freedom, 
or p{> Ax'^) = 0.05. Hence, while they are not definitively 
discrepant, there is still a suggestion of either increasing halo 



ellipticity, or increasing alignment of halo light and mass 
with luminosity. 

While one might be concerned about the model- 
dependence of these NFW and TIS predictions, we have 
verified that for red galaxies, changing the scale radius ra 
by 20 per cent in either direction changes the central value 
of fh by roughly 0.5(t, so the conclusions we derive from 
these results are not too strongly sensitive to errors in our 
derivations of r^. (For NFW profiles, that change in cor- 
responds to either a 20 per cent change in c(M, z) for a given 
halo mass, or a 25 per cent change in the halo mass using 
our assumed c(M, z) relation. For red galaxies, the central 
halo mass is generally known to this precision or better, as 
in Table |K|) For blue galaxies, the halo mass is not as well- 
determined, but neither is f — fa, so once again the possible 
systematic error in fh due to uncertainty in is within la. 
Furthermore, since the results with NFW and TIS profiles 
agree at the la level, it seems that they are not too sensitive 
to the exact form of the density profile assumed, as long as 
it is one that is steeper at larger radii, with the scale radius 
associated with this change increasing with luminosity as we 
have used here. 

Next, we analyzed the results assuming a power-law 
model for AE, which is approximately true for the scales 
under consideration, doing a joint fit for the amplitude and 
slope of AE, and for fh with /AE — /45 AE. In our case, the 
amplitude and a are just nuisance parameters, with the im- 
portant quantity being the best-fit value of fh ■ fh is slightly 
correlated with a due to the fact that the predicted value of 
/ — /45 can be related to a parameter combination including 
fh and a. Best-fit values of fh are shown in table |S] a is 
typically ~ 0.8-0.9. The fits were done individually for each 
lens-source sample combination, and the results were aver- 
aged over lens or source samples assuming Gaussian errors. 
As shown in this table, the results are similar to before in 
that (1) the results for the r < 21 sample are lower than for 
the other samples, (2) for blue lenses, the results are neg- 
ative for L3-L5 but positive (with large error) for L6, and 
(3) for red lenses, the results seem to imply some increase of 
ellipticity or of alignment of light and mass ellipticities with 
luminosity. However, as shown here, the values of fh implied 
for red lenses are lower than from the NFW and TIS anal- 
ysis from the non-Gaussian / — /45 determination. This is 
to be expected from our analysis in H2.3I and since the non- 
power law profiles are more realistic, the results in Table |K| 
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Table 5. Expected values of averaged over radial ranges described in the text, and over source samples for each lens sample, using 
SIS, NFW and TIS profiles with relevant parameters shown here. Results arc also shown averaged over luminosity and colour. 



Lenses 


SIS profile 

fh 


Mi8o (10"/i-i 


NFW profile 

Mq) (/i-lkpc) 


fh 


TIS profile 
Ts (/i-lkpc) fh 


L3, blue 
L4, blue 
L5, blue 


-1.1 ± 1.1 

-1.5 ± 1.0 
-1.1 ± 1.1 


6.1 ±2.8 
8.9 ±3.4 

9±8 


16 
19 
20 


1 7+2.5 
l./_2 y 

-9 1+*-'^ 

y-^-5.6 


180 
180 
250 


-Z.0_2 1 
9 c:+2.8 


L6, blue 
All, blue 


4 0+5.4 
^•'-'-3.6 

-1.1 ± 0.6 


165 ± 246 


81 


o O + 15.0 
-1 4+^-^ 


400 


-0.5 ± 1.3 


L3, red 


-1.3 ±0.5 


5.1±2.5 


15 




200 


-■-■^-1.3 


L4, red 


0.06 ±0.26 


13.0 ±2.7 


23 


0.3 ±0.9 


250 


0.59 ± 0.62 


L5, red 


0.31 ±0.32 


67 ±8 


52 


0.40 ±0.57 


350 


0.47 ±0.75 


L6, red 


1.2 ±0.5 


311 ± 82 


108 


1.7±0.7 


600 


2.7 ± 1.0 


All, red 


-0.06 ± 0.19 






0.60 ±0.38 




0.57 ±0.41 



Table 6. Results for from power-law profile fits described in the text for each combination of lens and source samples. All errors 
shown are Gaussian 68 per cent confidence intervals. 



Lenses 


r < 21 sources 


r > 21 sources 

fh 


LRGs 


all 


L3, blue 


-6.1 ±2.7 


-1.7± 2.0 


-3 ± 92 


-3.3 ± 1.6 


L4, blue 


-1.5 ± 1.8 


-0.33 ± 0.99 


0.0 ±4.2 


-0.6 ± 0.8 


L5, blue 


-0.29 ± 0.66 


-3.5 ±2.1 


-1.4± 2.0 


-0.7 ±0.6 


L6, blue 


13 ±9 


10± 11 


-0.4 ± 2.0 


0.5 ± 1.9 


All blue 


-0.7 ±0.6 


-1.0 ±0.8 


-0.8 ± 1.3 


-0.8 ±0.4 


L3, red 
L4, red 
L5, red 
L6, red 
All red 


-0.26 ±0.37 
-0.16 ±0.22 
-0.23 ± 0.22 
0.01 ±0.36 
-0.17 ±0.13 


-0.8 ±0.3 
0.11 ±0.18 
0.07 ±0.17 
0.21 ±0.08 
0.13 ±0.07 


-1.0 ± 1.0 
0.37 ±0.36 
1.7 ± 0.7 
0.88 ±0.42 
0.63 ±0.25 


-0.61 ±0.23 
0.05 ±0.13 
0.03 ±0.13 
0.21 ± 0.08 
0.10 ±0.06 



are probably more trustworthy. We note that in comparison 
with the SIS model results in Tabled the results with this 
power-law fit method averaged over all sources in the last 
column of TableEl(l) are diflerent, and (2) have significantly 
smaller errors. Considering that the power-law fit model has 
an additional free parameter (the power-law slope a) which 
should ostensibly lead to larger errors rather than smaller 
ones, and that the change from a. = 1 (SIS) to the best- 
fit values a = 0.8-0.9 cannot account for the difference in 
the results, these discrepancies require reconciliation. How- 
ever, it is important to keep in mind that the method used 
to determine fh for Table |S] correctly incorporates the non- 
Gaussianity of the measurement, whereas the power-law fits 
used for Table |S| have quoted errors from Gaussian propa- 
gation of errors on AE and /AE; consequently, the larger 
errors in Table|5]are more likely to be correct. Furthermore, 
since the use of the non-Gaussian error distributions weights 
the fh values in different radial ranges differently than the 
optimal weighting inherent in the c/ii^-minimization 

power-law fits, it is not surprising that the results from the 
two methods differ. 



5.6 Alignment with LSS 

Here we show results for the same tests using 800 /i~^kpc to 
2 /i~^Mpc in order to determine if there is any net alignment 



between the ellipticities of the hght and of local LSS (e.g., 
along filamentary structures). For these tests, we assume 
a power-law profile as is approximately observed on those 
scales, and use both the non-Gaussian determination and 
the power-law fit for fh that assumes Gaussianity of the 
errors. 

The non- Gaussian analysis, averaged over source and 
lens samples, yields a value of fh = —0.46 ± 0.34 for spirals, 
and —0.07 ± 0.09 for ellipticals. As mentioned previously, 
we expect that for spirals we should see some tendency to- 
wards negative fh on large scales, but unfortunately the er- 
rors are too large t o be able to confi rm this result at this 
time. For ellipticals, iKim et al] ll2002h showed that at least 
those that are the BCGs of clusters should have their ellip- 
ticities aligned with the cluster ellipticity, and therefore we 
may expect positive fh for these galaxies; this trend is not 
observed in this dataset. 

The Gaussian analysis (a power-law fit to the signal for 
amplitude and slope as nuisance parameters in addition to 
fh) on these scales yielded fh ~ — 0.42±0.28 for spirals, and 
—0.01 ± 0.05 for ellipticals. Hence, both analysis methods 
yield results that are similar, with a negative result that is 
not statistically significant for spirals, and a null result that 
places a constraint on the alignment of hght with LSS on 
0.8-2 /i~^Mpc scales for ellipticals. 
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6 DISCUSSION 

We can draw a number of interesting conclusions from the 
results in ^ First, several of our systematics tests have in- 
teresting implications. It is clear that in SDSS data, there 
are several problems that lead to difficulties in the measure- 
ment. The significant systematic shear, which seems to be 
related to the modeling of the PSF or the PSF correction 
itself, requires us to undertake additional complications in 
our analysis (the subtraction of / — /45 instead of comparing 
/ directly with the models). Second, small-scale photomet- 
ric problems lead to loss of sources preferentially along the 
lens major axis. 

Some of the systematics themselves have interesting sci- 
entific implications. For example, we have a clear detection 
of anisotropy of the distributions of satellites around lens 
galaxies, at least for the brightest red galaxies (—21 > Mr > 
—23). Since many of these galaxies reside in clusters, this 
may refiect that the BCGs (i.e., our lenses) arc aligned with 
the major axis of clusters, as foun d by Binggeli (1982) with 
44 AbeU clusters out to z = 0.1. [Fuller et al.f (ll999^ with 
poor, low redshift clusters. IWest fc Blakesleel f200(f) in 3 di- 
mensions for the Virgo cluster, and Kim ct al. (2002) using 
~ 300 clusters in SDSS data over a large range in redshift, 
0.04 < z < 0.5. This effect may be explained in terms of 
anisotropic infall into the cluster potential well along fila- 
ments. An upcoming work will explore this finding in more 
detail. 

Another interesting systematic is the anisotropy of 
sources due to magnification bias, of which hints were seen 
in the results, but which was not the focus of this work; 
future works with higher S/N on the anisotropy of lensing 
signal and magnification bias, and lower contamination by 
physically associated sources, should explore the relation- 
ship between the two in more detail, since they must be 
related for consistency to the average dark matter halo pro- 
files and alignment of light and mass ellipticities, and may 
provide another measure of fh that will allow for a more 
precise measurement of halo ellipticity. 

Our results for halo ellipticity are not conclusive. For 
spirals, we appear to have signs of anti-alignment of the el- 
lipticities of light and halo on a l-2a level (when averaged 
over luminosity) depending on the model used; for ellipti- 
cals, we see signs of a progression, with eh/e^ increasing 
with luminosity. If we trust the many theoretical predictions 
of dark matter halo ellipticity (which have been verified by 
many different types of simulations, both N-body and hy- 
drodynamic), then we must conclude that the trend we see 
refiects a change in the average alignment of the halo and 
light ellipticities with luminosity. 

We emphasize results with the non-Gaussian averaging 
over radial ranges rather than the power-law fits that ignore 
non-Gaussianity in the error distributions, and we empha- 
size those results using the NFW and TIS models rather 
than the SIS or other power-law models for (/ — /45)//h. 
Our preference stems from the fa ct that, while our res ults in 
iMandelbaum et all (l2005b^ and ISeliak et all (l2005all show 
that AE may appear to be a power-law on all scales shown 
there, 20 h'^kpc to 2 /i~^Mpc, according to the halo model 
that is used for the fits to the signal, the signal due to the 
central halo itself is fit more accurately with an NFW profile 
(or other truncated profile that falls below the power-law on 



several hundred /i~^kpc and larger scales), with the balance 
of the signal on those larger scales coming from satellites 
within the same halo for those galaxies that are in a group 
or cluster. If one assumes that the central halo itself has 
a power-law density profile, then tries to add in the signal 
from other galaxies in the group/cluster, the predicted sig- 
nal becomes too large to fit the data on large scales. Because 
of the change of the density profile with radius in these more 
realistic NFW/TIS models, the predicted signal anisotropy 
on the largest scales probed in this paper (~ 300 h^^kpc) 
is significantly smaller than the predicted anisotropy for a 
power- law profile (with larger error), leading to a signifi- 
cantly larger measured value of fh using these profiles. We 
have shown that our results are not very sensitive to even a 
20 per cent change in the measured r^, and are not too dif- 
ferent for the NFW or TIS, so they are not too dependent on 
the parameters of the density profile as long as that density 
profile roughly matches that of the expected one according 
to the halo model fits, with a scale-dependent power-law 
exponent. 

There is only one previous measurement of halo elliptic- 
ity using weak lensing. lHoekstra et alj (12004"). In that paper, 
the result was that fh = 0.77lg'2i. Our results, while more 
complicated due to the apparent luminosity dependence, do 
not seem at first glance to agree with this result. However, 
we must first ensure that we are comparing the appropri- 
ate things, since these two papers use different data and 
methodology. 

There are several differences in the data used for these 
papers. The SDSS data has redshifts, or at least photo- 
metric redshift s with (j{z) ^ 0.04, for all the lenses used; 
iHoekstra etall (|20o3) uses data from RCS, with redshift 
distributions only for both lenses and sources. Furthermore, 
since their lenses have 19 < ii < 21, and mean redshift of 
0.35, they correspond mainly to our L4-L6. In this work, 
we use colour infor mation to approximate ly separately el- 
lipticals and spirals; IHoekstra et al.l (|20o3) lack that capa- 
bility, so their results are for a mixture of the two. One 
might expect that if ellipticals have light and mass ellip- 
ticities ali gned , and spirals have them anti-aligned, then 
iHoekstra et al.l l|2o3) should have seen a null result for the 
mixture, but this is not the case. Indeed, if we average our 
results (including non-Gaussianity) from Tablej^over lumi- 
nosities L4-L6 and over colour, we find fh = 0.02 ± 0.21 
for the SIS because it does not take into account predicted 
radial-dependence of the results, or fh = 0.73 ± 0.39 for the 
NFW profile and 0.86 ± 0.43 for the TIS (note that these 
values are larger than the results averaged over L3-L6 be- 
cause L3 gives a significantly negative value, and they are 
positive because the signal for red galaxies entirely domi- 
nates over that of the blue galaxies) . The la tter two results 
do no t actually contradict the result from IHoekstra et alJ 
('2004*) , so even though we cannot say they constitute a 
definitive measurement of halo ellipticities, we do not have 
to worry about a significant discrepancy either. However, 
we note that this procedure of averaging the fh values over 
colour and luminosity, taking into account non-Gaussianity 
of error distributi ons, may not give th e same answer as the 
method used in .Hoekstra et aP ll2004l) . for which the signals 
themselves were all averaged and the maximum likelihood 
method was applied to the average signal. Furthermore, the 
maximum likelihood method itself is susceptible to bias due 
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to systematic effects such as the intrinsic ellipticity-density 
alignment that was observed to be present at some level 
(but taken into accou nt by a correction factor) in this work. 
iHoekstra et alj ll2004f) has a mean lens redshift of 0.35 and 
mean source redshift of 0.53, whereas this work has a lower 
mean lens redshift for all luminosities (though the bright- 
est bin is close to that of iHoekstra et al]|2004^ and lower 
mean source redshift for all but the LRG sample. Therefore, 
if there is some evolution in galaxy ellipticity with redshift, 
which has bee n proposed for cluste r ellipticities based on 
simulations in iHopkins et alJ i2005ri , with higher redshift 
clusters having higher ellipticity, then one might have ex- 
pected a larger detection in Hockstra ct al. ( 2004) than here. 

The two papers also fea ture different analysis methods. 
As argued in an appendix of lHoekstra et al] i2004) , they do 
not expect large systematic shear, though no measurement 
of /45 was presented to show that this is the case. However, 
systematic shear can only cause an und erestimate of /, not 
an overestimate. IHoekstra et alJ i2004) computed the the- 
oretical shear for a TIS model, then introduced anisotropy 
into the shear (i.e., not into k and then into the shear via 
derivation of the full 2-dimensional potential as done here) 
assuming that the anisotropy is equal to that in a SIS.* Since 
we have seen that the anisotropy for TIS approaches that of 
a SIS on small scales, but decreases with scale, the deriva- 
tion of fh = O.77I0 21 '® actually an under-estimate, since it 
assumes the SIS prediction even on larger scales. The fact 
that the SIS prediction for f / fh was used means that rather 
than compare our NFW/TIS predictions against theirs, we 
must compar e our SIS predic t ion of fh = 0.02 ±0.21 against 
the result of lHoekstra et al.1 (|20^, giving a 2.5cr discrep- 
ancy i f we take the errors at face value. The IHoekstra et alJ 
(|200^ paper derives errors on fh from the maximum like- 
lihood method, using the Ax'^ values. While this method 
typically gives reasonable confidence intervals even for non- 
Gaussian error distributions, the different error derivation 
between the two results nonetheless means that the appar- 
ent 2.5(T discrepancy should not be considered a cause for 
alarm. 

It is clear from this discussion that there are numer- 
ous differences between the two works, making the cause 
for the discrepancy uncertain, the primary contenders be- 
ing systematic contamination and physical differences of the 
samples used in the analysis. It is, however, also clear that 
we do not find a statistically significant alignment between 
light and dark matter. This may be expected from a the- 
oretical perspective, since th e presence of baryons te nds to 
round the dark matter halos iKazantzidis et al]|2004) . 

Future work to measure halo ellipticity with galaxy- 
galaxy weak lensing in the SDSS should have several fo- 
cuses. First, with time there will be more data, allowing us 
to increase the signal to noise of these measurements. Sec- 
ond, once there is more data, we will have the ability to 
characterize the systematics, such as systematic shear and 
the contamination by intrinsic alignments, better. Finally, 
we hope in the near future to have a catalog of clusters with 
ellipticities. Since we are finding hints of a positive detection 
of halo ellipticity in L6, the most massive lens sample, we 
hope for more definitive constraints using clusters as lenses, 

^ H. Hoekstra, e-mail communication. 



since they are more massive and the signal-to-noise on their 
lensing signal is higher than for galaxies. 

In conclusion, while we have not made a definitive de- 
tection of dark matter halo ellipticity using SDSS data, this 
work has made several contributions to future efforts in this 
field. We have identified some of the contaminants to such 
a measurement that may exist in all datasets, and there- 
fore should be considered in future analyses in this field. 
We have found suggestions of an anti-alignment of light el- 
lipticities with halo ellipticity on the 1-2g level for spiral 
galaxies (averaged over luminosity), and a suggestion that 
the light and mass ellipticity alignment in ellipticals is a 
function of galaxy luminosity, ranging from consistent with 
zero up to a 2g alignment in L6. Finally, we have attempted 
to find correlations between the ellipticities of light and LSS 
on 0.8-2 ft^^Mpc scales, and have placed constraints on such 
a correlation for ellipticals; results for spirals had fairly large 
errors so, while the results were negative in accordance with 
predictions, there is no clear detection. Future efforts with 
more data should allow us to make more definitive state- 
ments about these effects. 
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